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Abstract 

An analytic expression is proposed for the three-point function of the exponential 
fields in the Liouville field theory on a sphere. In the classical limit it coincides with what 
the classical Liouville theory predicts. Using this function as the structure constant of the 
operator algebra we construct the four-point function of the exponential fields and verify 
numerically that it satisfies the conformal bootstrap equations, i.e., that the operator 
algebra thus defined is associative. We consider also the Liouville refiection amplitude 
which follows explicitly from the structure constants. 



1. Introduction 

Since early 80's when the two-dimensional Liouville field theory (LFT) was recognized 
[1] as the effective field theory of the 2D quantum gravity considerable efforts has been 
directed at this area, especially for its relation to the string theory [1] (see e.g. refs.[2-4]). 
However despite significant progress in understanding the situation up to now the solution 
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to LFT is generally lacking. E.g., the structure of the operator algebra and the correlation 
functions in the general case are still unknown. 

Interest in LFT was renewed recently with the development of the alternative matrix 
model approach to the 2D gravity [5,6]. As the result of this new activity it was shown 
that LFT is able to reproduce some of the predictions of the matrix model approach, in 
particular the scaling behavior [7-9], the genus one partition functions [10] and some of the 
integrated correlation functions [11-15]. We would like to emphasize refs.[12] and [14,15] 
where the analytic continuation in the order of perturbation theory has been used. This 
approach is conceptually close to the way we make our guess about the structure constants. 

The content below is arranged as follows. In sect. 2 the LFT on a sphere is introduced 
and some notations are defined. In sect. 3 we propose an exact expression for the LFT 
three-point function and discuss some of its properties. The Liouville reflection amplitude 
is also introduced here. It should be stressed that the argvmicnts of this section have 
nothing to do with a derivation. These are rather some motivations and we consider the 
expression proposed as a guess which we try to support in the subsequent sections. This 
guess appears quite natural^ and might even be thought obvious to those concerned with 
the problem. What we believe does promote it to a step forward are the various tests 
performed in sects. 4-7. 

Sect. 4 contains some calculations in the semiclassical limit. The physical content of 
the Liouville reflection amplitude is discussed in sect. 5 and in sect. 6 we make a check of 
this amplitude by means of the thermodynamic Bethe ansatz for the high-temperature 
sinh- Gordon model. 

The exact LFT structure constants together with an effective numerical algorithm for 
the conformal block permit us to compute numerically the four-point function of the LFT 
exponential flelds. We perform this in sect. 7 and verify that the four-point function satisfles 
the conformal bootstrap equations, i.e. the necessary conditions of the associativity of the 
operator algebra. Sect. 8 contains some considerations about the classical Liouville action 
and the related problem of accessory parameters. 

2. Liouville field theory 

Local properties of the LFT are derived from the Lagrangian density^ 

C=^{dacl>f + fie''^ (2.1) 



^ Compare (3.14) below with the formulae for the structure constants of "minimal 
CFT"[16] and other "rational CFT"[17] which were taken as "architectural prototypes" 
for (3.14). 

^ It is conventional to add also another term ^R\/g (j) to the Liouville Lagrangian 

density, with g and R being arbitrary "background metric" and associated curvature, and 
with the parameter Q adjusted to ensure that all physical quantities be independent of a 
particular choice of this "background" (see [8,9]). However, it is always possible to choose 
a speciflc background which is flat everywhere except for few selected points; the above 
term translates then into appropriate "boundary" terms, as in (2.7) and (2.33) below. 



2 



where b is the dimensionless Liouville couphng constant and the scale parameter /j, is usually 
called the cosmological constant. Below we often use the complex euclidean coordinates 
z — xi + ix2; z = xi — ix2 and denote d = d/dz; d = d/dz. The Liouville field (j){z, z) is 
not exactly a scalar but varies under holomorphic coordinate transformations z ^ w{z) as 



(p{w,w) = (p{z,z) 



Q 



log 



dw 



dz 



where 



Q = b+l/b 
The holomorphic Liouville stress tensor 

T{z) = -{d(Pf + Qd^<p 
T{z)^-{d4>f + Qd^4> 

ensures local conformal invariance [18] of LFT with the Liouville central charge 

CL = l + 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



To define LFT globally one has to specify boundary conditions. The LFT on a sphere 
corresponds to defined on the whole complex plane with the following asymptotic be- 
havior at \z\ oo 

(j){z,z) ^ -Q\og{zz) + 0{l) at 1^1 ^oo (2.6) 

For specific calculations it is sometimes convenient to set this asymptotic behavior by 
considering LFT on a large disk F of radius i? — > oo and adding a boundary term to the 
Liouville action 

Al = ^ I [{da(l>f + A^T^ie^^n d^x+-% f (f)dl + 2Q^ log R (2.7) 
47r Jr TTit Jgr 

The last constant term is introduced to make the action finite at — oo. This type of 
boundary condition is conventionally called the background charge —Q at infinity. 

Exponential Liouville operators 

Va{x) = e2"^(^) (2.8) 

are the spinless primary conformal fields of dimensions 

Aa = a{Q - a) (2.9) 

Note that the field Vq_q, has the same dimension as (2.8). These two fields are closely 
related and we shall call Vq_q, the refiection image of and vice versa. In particular the 
perturbation operator Vb in (2.1) have dimension 1 together with its refiection image V1/5. 
The case a = (5/2 is degenerate. Here we have two primary fields 

^Q/2(a;) = eQ'^(") (2.10) 
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and 

Uq/2{x) = V^{x)\^^Q/, = 0(a;)e^^(-) (2.11) 

of dimension Q^/4:. The last field (2.11) is called sometimes the puncture operator in LFT 
[11]. 

The n-point function of the exponential fields on a sphere is formally defined as a 
functional integral 

Ga„...,aAxi,-..,Xn) = j K., (xi) . . . (x„)e-^^['^]D0 (2.12) 

over the fields 0(x) with the boundary condition (2.6). Note that in the definition (2.12) 
we do not divide the right hand side by the zero-point function Zq. 

The scale (//) dependence of any correlation function (2.12) [7-9] 

Qa„...,aSxu- • • , ^n) = (Tr/x)^^-^ (o^i, . . . , a;„) (2.13) 

(with -Fai,...,a„ (3:^1, • • • T^n) independent on jj) is easily derived from the action (2.7) and 
the operator product expansion 

(i){z,z)V^{x) = -alog\z-x\^Va{x) + ... (2.14) 

Sometimes it is convenient to consider the n-point functions Qkl^,...,a„{xi, ■ ■ ■ ,Xn) with 
fixed area 

A= f e^'"^(fx (2.15) 



These observables are related to (2.12) as 

ga,,...,aAxi,...,Xn) = J. Jo^i, . . . , a;„)e-'^^— (2.16) 

so that 

^.„...,a„ (XI, . . . , X J - ^ ^ J F ((g - E «^) /b) ^ ' 

Conformal invariance restricts to some extent the x-dependence of the correlation 
functions [18]. In particular the three-point function is specified up to an ^-independent 
constant C{ai,a2,as) 

^ai,a2,a3 (^1)^2,3^3) = 1^^231^'^' \x3l\^'^^ C{ai, a2, a^) (2.18) 

where 71 = A^, - - A^g, 72 = A^^ - A^g - A^i, 73 = A^g - A^i - A^^ and here 
and below we denote The four-point function can be reduced to a function 

of only one coordinate variable, the projective invariant of the four points xi, . . . ,a;4 

a;=^i^ (2.19) 

4 



as follows 

^ai,...,Q;4 (-^l ) • • • ) •^4) — _|_ _|_ I "^ij I ^ai ,02,0:3,0:4 ("^J (2.20) 



i<j 



where 712 = 713 = 0, 714 = -2Aai, 724 = + A^^ - A^^ - A^^, 734 = A^i + A^^ - 
Aas - Aa4 and 723 = A^^ - A^i - Aa^ - A^g. Function ^01,02,03,04(3^7^) satisfies the 
following symmetries (sometimes called the crossing symmetry relations) [18] 



^01,02,03,04(^5^) — ^01,03,02,04(1 x^l x) — \x\ "1(^0,^^0.^^0.3^0,2(1/3^,1/3^) (2.21) 

In principle one can separate two pairs of operators in the four-point function, say 
Vai{xi)Voi2{x2) and Vo3(a;3)yQ,4(3^4) and represent ^ai,...,a4(a;i, • • • , 2:4) as a sum over in- 
termediate physical states [18]. As was established long ago [2,3] the physical LFT space 
of states A consists of a continuum variety of primary states corresponding to operators 
Va with 

a=^ + iP (2.22) 

(P is real) and the conformal descendants of these states. The corresponding expression 
for the four-point function looks as 

,02,03,04 ('^J '^) Q 

/ C{ai, a2, Q/2 + iP)C{a3, a^, Q/2 - iP)\J^{A^^, A, x)\^dP 
—00 

(2.23) 

Here C{ai,a2,as) are the structure constants of eq.(2.18) while J-'{Acn, A, x) is the so- 
called conformal block [18] which sums up all the intermediate descendant states of a given 
primary one. The conformal block is determined completely by the conformal symmetry 
of the theory. It depends on the corresponding central charge (2.5), on the dimension 



n2 

A = + (2.24) 

of the intermediate primary state and also on the "external" dimensions Aq., i= 1, . . . , 4. 
Unfortunately up to now this function is not known in a closed form. However it is 
straightforward to evaluate it as a power series in x [18]. Also there are very efficient 
algorithms for its numerical computation [19]. 

One important remark about Eq.(2.23) is in order. As we will see (sect. 3) the structure 
constant C{ai, 02, ct) exhibits infinitely many poles in the variable a, their positions being 
dependent on ai,a2- As one changes ai,a2 and q;3,q;4 in (2.23) the associated poles of 
the integrand move around and some of them can cross the integration contour QmP = 0. 
Therefore the decomposition (2.23) can be taken literally only if ai,a2 and 0:3, 0:4 take 
their values within certain domains, namely 

|ai - Q!2| < Q/2; \Q - ai- a2\ < Q/2; 
\as — a4\ < Q/2; \Q — — < Q/2. 
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(here we consider the case of real a's). Otherwise the four-point function (2.23) has to be 
understood as analytic continuation away from the domain (2.25). In the course of this 
continuation some poles of C(ai, 02, (5/2 + iP) and/or of C{a^, 04, Q/2 + iP) cross the 
real P axis and the r.h.s. of (2.23) acquires additional "discrete" terms associated with 
the residues of the integrand at these poles (this seems to agree with qualitative analysis 
in [11]). We will discuss this phenomenon in greater details elsewhere. 

Of course one can pick up another partition of the four operators into pairs and 
obtain an expression similar to (2.23) summing up over the intermediate states in the 
corresponding channel. The resulting four-point function must be the same. In other words 
the four-point function defined as in eq.(2.23) has to satisfy the crossing symmetry relations 
(2.21). From this point of view these relations appear as a set of non-trivial restrictions 
for the structure constants C(ai,a2,Q;3) which express the associativity of the operator 
algebra and are known as the conformal bootstrap equations [18]. We discuss more about 
this point in sect. 7. Representations like (2.23) can be written down for higher multipoint 
functions as well. However they involve the multipoint conformal blocks which are much 
more complicated objects then the four-point one. Fortunately, the corresponding crossing 
symmetry relations are not expected to bring up any new restrictions on the structure of 
the operator algebra [18]. 

Let us finish this section with few words about the classical limit 6 ^ of LFT. Here 
it is more convenient to use the field 

^ = 26(/) (2.26) 

which becomes a classical Liouville field in this limit. Its dynamics is governed by the 
classical action 

'5Liouv[¥'] = fc'^L[0] at 6^0 (2.27) 



Explicitly 



d^x + c^oo + 21ogi? (2.28) 



where 

^00 = 77^ / Vdl (2.29) 

The field (p{x) satisfies the classical Liouville equation 

dd^ = 27r/x6^e'^ (2.30) 

and locally describes a surface of constant negative curvature — 87r//6^. 

The leading (exponential) asymptotic of the n-point function (2.12) in the classical 

limit 

Oa„...,aA^l, ...,Xn)-^ exp (^'^S^'^'^^ (2.31) 

is governed by the classical Liouville action on an appropriate solution to the Liouville 
equation. Note that the insertion of any Vq,. (xi) in (2.12) affects the classical field dynamics 
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only if the corresponding ctj is "heavy", i.e. if cuj = r]i/h is of order 0{h~^). Technically 
one has to distinguish two possibilities. If 77^ > 1 a classical solution with negative 
curvature exists and in (2.30) 

-^if.L.rj^Ca^l^ ■■■,Xn)= /SlIouv {x\xi, Xn)] (2.32) 

where <^^^^...^^^(a;|a;i, . . . , Xn) is a solution to (2.30) with the following boundary conditions 

V{z,z) = -2log\z\'^ + 0{1) at |2;|^oo 

(p{z, z) = —2r]ilog\z — Xi\^ + 0{1) at z ^ Xi 

This field configuration is singular at z Xi. Therefore it is better to cut out small disks of 
radius around each point Xi and define the regularized Liouville action on the remaining 
part r of the complex plane 

SuouA<f] = ^j^ lida'pf + STTfib'^e'^ d^x+<foo+2logR-Y,{vm + '^Vi^ogei) (2.34) 
Here <^oo is as in eq.(2.29) while the boundary terms with 

<Pi = 7^ f Vdl (2.35) 

arc added for each small circle to ensure the behavior (2.33) near Xi. Also we include some 
field independent terms such that Suonv is finite and independent on at ej —>■ 0. 

At ^ ?7i < 1 there is no solution to (2.30) and (2.33) with negative curvature. In 
this case it is relevant to consider the n-point function (2.17) of fixed area. The leading 
classical behavior is again governed by the classical action 

€t..,aM^ • • • '^n) - exp (-^S^f) (2.36) 
but now one has to impose the constant area constraint 

e'^cfx (2.37) 

r 

and solve the positive curvature Liouville equation 

dd^='^(j2vr-l)e^ (2.38) 



with the same boundary conditions (2.33). The corresponding constant area Liouville 
action is again (2.34) but without the cosmological term /ib'^A so that 

-^Liouvb] = S^f [ifi] + i^b^ e^£x (2.39) 



Note that the same fixed area calculation can be performed in the first case J^r/j > 1 
too. The integral over A in eq.(2.16) is now dominated by the stationary point A = 
{YlVi ~ which corresponds to the classical solution of eq.(2.30). 

Suppose now that the multipoint function contains several "heavy" operators with 
cti — Vi/b cind also a number of "light" ones, i.e. the Liouville exponentials with aj — ajb 
of the order 0{b). In the classical limit these "light" fields influence neither the classical 
solution nor the one-loop correction. Therefore one expects that as 6 — > 

Saib,...,aib,m/b,...,r]„/b{yi, ■ ■ ■ , Vl, Xl, . . . , Xn) ^ j-j- ^aj^ni,...,r,i{yj\xi,...,Xr,) (2.40) 
Orji/b,...,r)r,/b{xij ■ ■ ■ j^n) 

where (Pni,...,rii {x\xi, . . . , Xn) is the classical solution for the "heavy" operator configuration. 
Special care is required if the number of the "heavy" operators is less then 3. In this case 
the functional integral (2.12) has zero modes which have to be integrated out explicitly 
even in the classical limit. In particular, if there are no "heavy" operators at all the relevant 
fixed area classical solution is the 2D sphere metric of area A 

ifo{z,z) = log ^ (2.41) 
7r(l + zz)^ 

This solution has to be integrated over its SL{2,C) orbit parameterized by four complex 
numbers a, b, c, d with the constraint ad — be = 1 

ipQ{z,z\a,b,c,d) = log— ^ ,„ , — — (2.42) 

TT{\az + op + \cz + d\^y 

This leads to the following expression for the n-point function of "light" operators in the 
classical limit [20] 

g^f6,....a„6(^i.---.^n) e'^^^-^^^M^^'^cA)^^^^^ d) (2.43) 

^0 i=l 

where Z^^^ is the fixed area zero-point function (the partition function of the sphere) and 
(i/x(a, 6, c, d) stands for the invariant measure on SL{2, C), 

dii(a, b, c, d) = Ad^a d% d^c d^d 5^^^ (ad -bc-1). (2.44) 



3. Three-point function 

Naively one can try to expand the A'^-point function (2.12) into a perturbative series 
in the cosmological constant 

oo 

Qoci,...,aN{xi,...,XN) = X]^S-:«iv(^l'---'^Af) (3.1) 

n=0 



where 



6?£La. = / {VaA^i)...V^^{xNmu^)...V,{Un))d'u^...cfun (3.2) 

and (. . .) denotes the functional integral over a free field (i.e., (2.12) at = 0) 

= n - 



(3.3) 



i>j 



However it is well known that in LFT expression (3.1) in general does not work. The 
problem is that the free field functional integral (3.3) for the n-th term (3.2) matches the 
spherical boundary condition (2.6) only if 



TV 



^ai = Q-nb (3.4) 



Following refs. [12-14] we shell interpret eq.(3.4) as a kind of "on-mass- shell" condition. 
Namely, Qai,...,aN (^i; ■ ■ ■ ^^n) exhibit a pole in the variable a = ^ ctj every time eq.(3.4) 
is satisfied forn = 0, 1, 2, . . ., the residue being specified by the corresponding perturbative 
integral 

res Qai,...,aMixi,---,XN) = Q^o^,...,aMi^^^---^^N) ^ ^ ^ (3.5) 



ai=Q—nb 



ai=Q—nb 



It seems unlikely that the on-mass-shell condition (3.5) alone is enough to determine 
the iV-point function. The analysis is hampered by the quite complicated and in general 
unknown analytic structure of the multipoint perturbative integrals (3.2). The situation is 
simplified in the three-point case N = 3. The a;-dependence of the on-mass-shell integrals 
turns out to be the same as that of the conformal three-point function (2.18) 

^Sa2,a3(^i'^2,a;3) „ = |a;i2|^^^ ^23!^^' ksil^^^ ^(ai, "2, "3) (3.6) 

y ^ on=Q—nb 

with the same 71, 72 and 73 as in eq.(2.18), while /n(Q;i, 0:2, cus) has been carried out 
explicitly in ref.[21] 

, ( -^^Y n?^i li-jh') 

ln[Oti,a2,OL2) — —, — TTTr ^ (3.7) 

Y[lZl[i{2arh + kh^)^{2a2h + kh^)^{2a^h + kh^)] 
Here and below the standard notation 

7(0;) = r(x)/r(i - x) (3.8) 
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is used. The on-mass-shell condition now reads 



res C{ai,a2,a3) = In{ai,a2,a3) 

ai=Q—nb 



(3.9) 



At this step we need to introduce some special function T(6, a;). Below we consider b 
as a parameter and suppress it in the notation. In the strip < Hex < Q function T(a;) is 
defined by the integral representation 



X 



sinh' 



X 1 ^ 



sinh ^ sinh 

2 26 J 



(3.10) 



From the definition it is clear that 

T(a;) = T(Q - x) 
T(g/2) = 1 

and that T(a;) is "self-dual", i.e. remains unchanged if b 
functional relations are derived 



(3.11) 

1/6. Also the following 



(3.12a) 
(3.126) 



T(a; + 6) =7(6a:)6^-2^^T(a;) 
T(x + 1/6) = 7(x/6)62^/^-^T(x) 

Using (3.10) and (3.12) it is easy to verify that T{x) is an entire function of x with zeroes 
located at a; = —m/b — nb and x = (m -|- l)/6 -|- (n -|- 1)6, where m and n run over all 
non-negative integers. Below we also use the notation 

dT{x) 



dx 



(3.13) 



x=0 



With the relations (3.12) it is straightforward to verify that 



C{ai,a2,as)= 77/^7(6^)6^ '^^ 



iQ-^ai)/b 



X 



ToT(2ai)T(2a2)T(2a3) 



T(ai + a^ + as — Q)T(ai + 0:2 — a3)T(a2 + — ai)T(as + ai — 02) 

(3.14) 

satisfies the on-mass-shell condition (3.9). We propose this expression as the exact three- 
point function in LFT. 

As a function of a = l^i=i "^i expression (3.14) has more poles then predicted by 
eq.(3.9). They appear at o; = Q — m/b — nb and at a = 2Q + m/b + nb for any pair of non- 
negative integers m and n. It seems suggestive to note that the corresponding residues are 
related to the more general multiple integrals, also evaluated explicitly in ref.[21]. Namely 

res Qai,a2,a3{xi,X2,X3) = 

cti=Q — m/b—nb 



mini 



// 3 m n \ 

( n (^^) n ^i/bM n ^^ivj) ) d^u^ . . . d^u^d^v^ ...d^ 
\i=i I 



(3.15) 
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where the "dual" cosmological constant p, is related to /x as 

nfi^il/b') = {n^i^ib')^ (3.16) 

The whole expression (3.14) is self-dual in the sense that it is invariant under the substi- 
tution 6 — > 1/6, |U ^ /i. 

The three-point function (3.14) reveals some interesting features. In particular the 
reflection a — Q — ct of each of the three operators introduces the Liouville reflection 
amplitude S{P) 

C{Q — CKi, 0:2, 0:3) = C{ai, a2, a3)S{iai — iQ/2) (3-17) 
which reads explicitly 

a( m _ (^,.^(h^)\ -2^^/^ ^(l + 2^P/5)^(l + 2^P6) 

S{P) - - {n^,^{b )) r(l - 2iP/b)T{l - 2iPb) ^^-^^^ 

This function is discussed more in sect. 5. We shall also use the notation 

G{a) = S{-^a + ^Q/2) = ^ (2 - 2a/b + 1/b^) ^'-''^ 

associating G{a) with the two-point function of operators Va{xi)Va{x2) in LFT. 

There are several different classical limits of the three-point function (3.14) dependent 
on how the parameters behave as 6 — 0. We consider here two cases. First, let all the 
three operators be "heavy", i.e. ai = rji/b aX b ^ with rji fixed. The leading asymptotic 
is 

C(cti, a2, cts) ~ exp 

,''12, Vs) (3.20) 

where 

S^^^^ ivi , = CZ^r]i- 1) log(7r/x6^) + F(r/i + r/s + r/s - 1) + F{r]i + r]2 - V3)+ 

F{m + ^3 - ^i) + i^(% + ^1 - m) - ^(0) - F(2?7i) - F(2?72) - i"(2?73) 

(3.21) 

and we have denoted 

F{r]) = / \og-i{x)dx (3.22) 

Jl/2 

Note that S^^^^rji, 772, 773) vanishes if Xli=i Vi = 1- 

In the opposite case of three "light" exponentials with ojj — bai of the order 0(b) it 
is more relevant to consider the fixed area three-point function (2.17). For this we find at 

C^^\bauba2,ba,) = [-] x 

r(cri -I- 0-2 + tJs - l)r(cri + 0-2 - 0"3)r((72 + ffs - tTi)r((T3 -I- (71 - (72) 

r(2cTi)r(2a2)r(2(T3) 

(3.23) 
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where C is the Euler's constant. 

We present also the classical limit of the fixed area two-point function G^'^^ (a), related 
to (3.19) as in eq.(2.16). In the "heavy" limit a = rj/b, 6 — > the leading exponential 
asymptotic 

G(^)(V6)~exp(^-l5('=')(r/)) (3.24) 
is governed by the classical action 

S^^^^ {r}) = (1 - 2r}) (^log ^ + log(l - 2ri) - 1^ (3.25) 

The fixed area two-point function of two "light" operators with a — ba reads at 6 — > 

(^^) ^ ±\ (3.26) 

^ ' \n J V2^63(2cT-l) ^ ^ 



4. Classical limit 

The asymptotic (3.20), (3.21) has to be compared with the classical action (2.34) 
evaluated at the solution (^^^^^2,r/3(-2|^ci, a^2, a^s) of the classical Liouville equation (2.30) 
with three singular points xi, X2 and x^. The boundary conditions (2.33) at these points 
are automatically imposed by the boundary terms in (2.34). In the three-point case the 
solution </'?7i,?72,?73(-2|a;i, X2, X3) can be found explicitly in terms of hypergeometric functions 

Vm,V2,mi^\^^^^^^^3) = -2log[ai'il;i{z)'il;i{z) + a2'ip2{z)'ip2{z)] (4.1) 

Here 

M^) = iz- a^i)"! (z - x^y-^'-^'iz - xzY^2Fi{m +V3-V2,Vi+V2+V3-^,'2Vi,x) 

^1;^{Z) = {Z - Xiy-'^'iz - X2)'''+'''~\z - XsY-'^'X 

2Fi(l + 772 - ?7i - ?73, 2 - r/1 - r/2 - 773, 2 - 2?7i, x) 

(4.2) 

and we have denoted 

X = (4.3) 

[z - a;2)a;3i 

while the constants ai and 02 read explicitly 



2 7r/x6^ 

"1 ^ 1^ |4r,3+4r?i-2 I |2-4r,2, \2-4r,i 
\Xl3\ \Xl2\ \X23\ 

72(2771)7(772 -hr/3-?7i) 
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(4.4) 



and 

l^isl^ (1 - 2r]iYai 



a. = - . (4.5) 



Now we have to evaluate the integral (2.34). Calculation of the Liouville action can 
be simplified by the following trick. Since <frji,...,'n„{^\^ij ■ • ■ j^n) is an extremum of (2.34) 
we have 

-K—Suouw[<Pm,-,Vu{^\^i^---^^n)] = -<yfi - 4r7i log ei (4.6) 
where (pi is defined in eq.(2.35). Near the singular point Xi the solution behaves as 

Vvi,-,vu{^\^u---,xn) = -2r]ilog\z - Xif + Xi + O {\z - Xif-'^'^') (4.7) 
(we suppose here that all ry^ < 1/2). Therefore in the limit — > 

d 

— 'S'uouv = (4.8) 

This equation implies that the form 

n 
i=l 

can be integrated, defining Suouy up to a constant independent on rji. To fix this inte- 
gration ambiguity we note that the Liouville action (2.34) can be explicitly evaluated if 

Er=i Vi = 1 

'S'Liouv|^^. = l = ^'^ViVj^Og\Xi- Xjl'^ (4.10) 
i<j 

In the three-point case the coefficients are easily derived from the explicit solution 
(4.1-5). E.g. 



Xi = -(1- 27/1 ) log 



— log TTIlb^ — 

"^23 (4.11) 
livi + V2 + V3- l)7(m + V2- Vahivi + V3- m) 

72(2771)7(772 + 773 - Til) 



while X2 and X^ are obtained from (4.11) by obvious permutations of 77^ and Xj. Integrating 
the form (4.9) we find 

'S'Liouv = C^'Hi - l)log7r/i6^ + {5i + 82- (^3)log|a;i2|^ + (^2 + (^3 - h)^og\x2z\^+ 
{h + 82) log |xi3p + F(77i + 772 + 773 - 1) + F(77i + 772 - 773)+ 

Fim +V3- Vi) + F{V3 + V1- m) - F{2r]i) - ^(2772) - ^(2773) - Fo{xi) 

(4.12) 
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where Si = r]i{l — rji), function F{rj) is defined by eq.(3.22) and Fo{xi) is the integration 
constant independent on 77^. Comparing with (4.10) we find FQ{xi) = F{0) in complete 
agreement with the asymptotic (3.21) of the proposed exact three-point function. 

Similar and even more simple calculations with the positive curvature Liouville equa- 
tion (2.37) support the classical limit (3.25) of the two-point function. 

One can also check that in the case of three "light" exponentials asymptotic form (3.23) 
of the proposed three-point function is in agreement with explicit semiclassical calculation 
through the equation (2.43). In the case n = 3 the x-dependence of the integral in (2.43) 
can be easily isolated with the help of its SL{2,C) transformation properties, so that this 
integral takes the form 

ki2|'"1a;23|'"Ha;3iP"^ I{(Ti,a2,a3), (4.13) 
where ui = ai — a2 - 0-3, 1^2 = cr2 — 0-3 — (7i, vs = os — ai — (72, and 

/((7l,(72,(73) = I — ) X / ((71, (72, (73); 



n 



/((7i,(72,(73) = J {\bf + \d\Y^'''{\a + bf+\c + d\Y^'''{\af + MY^"'dfi{a,b,c,d). 

^(4.14) 

It is this factor /((7i, (72, (73) which, being multiplied by the partition function Zq^"^ (see 
(2.43)), is to agree with the asymptotic formula (3.23). To evaluate the integral / in (4.14) 
it is convenient to use complex coordinates ^1,^2,^3 on the group manifold of SL{2,C) 



related to a,b,c,d as 



b ^ a + b ^ a , . -.^ ^ 

«' = 5' ^^ = 7^ = ? f^-^^' 



In this coordinates the invariant measure on SL(2, C) takes the well known form 
and the integral in (4.14) simplifies as 



(i//(a, 6, c, d) = — -^rxr;: T\vl^ (4-16) 



/(fTl, (72,(73) = j d%d%d% 



\Cl2\-'-"'^^\C23\-'-^'''\^3l\-^-^''H^ + I6|')-'"U1 + I6|')-'"^(l + 161')-'"^ 

(4.17) 

Now, it is straightforward to verify that this integral is invariant under the SU{2) subgroup 
of SL{2, C). Namely, the form of this integral does not change if one substitutes 

^ ■ ""^'^^ |a|2 + |6|2 = l. (4.18) 



-b$,i + a 

Using this symmetry one can set, say, ^3 = 00 and write the integral (4.17) as 

/((7l,(72,(73)= 7r j d2^id2^2|^,-^2|-2-2-3(l + |^^|2)-2a,(;L+|e2|')-''^^ (4.19) 
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This integral can be evaluated explicitly with the result 

~ 3r((Ji + (72 + CT3 - l)r(cri + (72 - Cr3)r(cr2 + CTs - CJi)r(cT3 + CTi - CT2) 

I{a,,a,,as)=^ 1X2^1)1X2^^)1^ ' 

(4.20) 

in agreement with (3.23). 

5. Reflection amplitude 

Consider LFT on an infinite flat cylinder of circumference 27r with the cartesian co- 
ordinates xi, X2 (as before z = xi + 1x2, z = xi — 1x2) and let us interpret the coordinate 
X2 along the cylinder as the (imaginary) time while xi ~ xi + 2?! be the space coordinate. 
The holomorphic Liouville stress tensor (2.4) allows one to construct two copies (right and 
left) of Virasoro algebra Vir and Vir with the central charge (2.5) 



[Lm, Ln] = (m - n)Lm+n + - rn)5m+n 

[Lm, Ln] = (m - n)Lm+n + ^^("^^ " 'rn)5m+n 

where the operators L„ and L„ appear in the expansion of the stress tensor 

CO 

n=— 00 

00 



(5.1) 



(5.2) 



and act in the space of states A of LFT on the constant time circle X2 = const. In particular 
the Hamiltonian 

H=-^+Lo + Lo (5.3) 

generates translations along the time X2. 

The space of states A is classified in the highest weight representations of Vir (8) Vir 

A = ®pAp (5.4) 

Each conformal class Ap contains a primary state vp which satisfies 

LnVp = LnVp — at n > 

Lovp = Uvp = {Qy4: + P^)vp ^^'^^ 

(so that the energy of fp is 2P^ — 1/12) and its descendants generated by the action of 
and Ln with n < on vp. Right and left generators Ln and L„ commute and therefore 
Ap has the structure of a direct product of right and left modules. 
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To get more of an idea about A take the "zero-mode" of the Liouville field 0(a;) 



^o = J%{x)^ (5.6) 

and consider the region (f)o — oo in the configuration space. Here one can neglect the 
exponential interaction term in the LFT action and consider (f){x) as a free massless field 
which can be expanded as usual in the free field oscillators 

cf,{x) = cf>o-V{z-z) + y f!^e-^ + !^e--^") (5.7) 

Here 

is the momentum conjugate to the zero-mode and 



r- - ] -—^ 



(5.9) 



(5.10) 



The Virasoro generators are represented at 0o — — oo as follows 

L„ = ^ akan-k + (2P + inQ)an n 

fc>0 

and the same for L„ with substituted by an- These operators act in the space of states 

Ao = £2{-oo < (j)o < oo) ^ (5.11) 

where J-' is the Fock space of the oscillators a^, an with n E Z, n 0. The space J-' consists 
of the Fock vacuum |0) annihilated by all a^, an with n > and the states generated by 
the action a„ with n < on |0). 

Let '^g[(l){xi)] be the wave functional of any state s e A. The 0o — —oo asymptotic 
of this wave functional is naturally associated with an element of (5.11). Prom (5.10) one 
observes that the zero-mode plane wave exp(2iP(^o) times the Fock vacuum |0) behaves 
under (5.11) as the primary state vp. It is clear that the correct wave functional [0(a;i)] 
of Vp contains at 0o —>■ —oo also a reflected wave exp(— 2iP(/)o) |0), i.e., 

[0(a;i)] ~ (e^^^-^^ + 5(P)e-2^^'^°) |0) at 0o ^ -oo (5.12) 
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with some reflection amplitude S{P) dependent on a more complicated dynamics in the 
region of 0o where the exponential interaction term is important (see [11]). More generally, 
if sp is some state of 

^sAH^i)] ~ + ^(P)e-'^^'^°) \s) , 00 ^ -OO (5.13) 

where |s) G and S{P) is now a unitary operator in J^. In particular S{P) is the eigenvalue 

of S{P) on the Fock vacuum |0). 

It is important to note that at a given P all the matrix elements of the operator S{P) 
are in fact determined by the conformal symmetry of LFT up to an overall multiplier. For 
example, applying L_i in the form (5.10) to the wave functional (5.12) we obtain 

O — HP 

SiP)a_, |0) = SiP)^^:^a_^ |0) (5.14) 

For more complicated Fock states the calculations are more involved but always reduce to 
linear algebra and permit us to restore S{P) uniquely up to S{P). 

The following features are readily established. First, is a tensor product of the right 
and left modules (generated by an and an respectively) and 

S{P)^S{P)sr{P)®sl{P) (5.15) 

where sr{P) and sl(-P) act independently in the right and left Fock spaces, sl{P) being 
isomorphic to Sl{P) under — * a^. Next, s^(P) commutes with Lq and therefore act 
invariantly at each level (i.e., at every eigenspace of Lq). Moreover, it can be argued that 
sr{P) commutes with the infinite series of the "quantum KdV integrals of motion" [22,23] 
built of the higher powers of the stress tensor. In other words, it has the same eigenvectors 
as the quantum KdV transfer matrix constructed in ref . [24] . 

However, generally Si?(-P) is not known in a closed form. We quote here its matrix 
elements at the second level spanned by the Fock states |1, 1) = a^.^ |0) and |2) = a_2 |0) 



D{P)sRiP) |1, 1) = (8P3 + (6Q2 _ 2)P + iQ{2Q^ + 1)) |1, 1) + 4iPQ |2) 
D{P)sr{P) |2) = -8iPQ |1, 1) + (-8P3 _ (6Q2 _ 2)p + iQ{2Q^ + 1)) |2) 



(5.16) 



where 

D{P) = (2P - iQ) (2P - i{b + 2/6)) (2P - i{2b + 1/6)) (5.17) 

The overall factor S{P) in (5.15) is not prescribed by the conformal symmetry and 
has to be recovered separately. It is easy to evaluate it in the semiclassical limit 6^0. 
Suppose that P is also small of the order of 0{b). Then even in the region of 0o where the 
interaction is significant one can neglect the oscillators in (5.7) and study only the dynamics 
of the zero- mode 0o- In this approximation, known as the minisuperspace approach [20] 
in LFT, the Hamiltonian (5.3) is substituted by 
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and the reflection amphtude (5.12) appears as 



-2iP/fe r (1 + 2iP/h) 

r (1 - 2iP/h) ' 







(5.19) 



Comparing this with (3.18) we find it natural to propose the expression (3.18) as the exact 
reflection amphtude S{P) in LFT. In the next section we verify this suggestion numericafly 
using the thermodynamic Bethe ansatz technique. 

6. Thermodynamic Bethe ansatz 

In this section we consider the sinh-Gordon model on a circle of circumference R with 
periodic boundary conditions. The problem is defined by the action 



^shG 



j dx2 J dxi 



1 2 

— {da(f>) + 2/1 cosh 60 



(6.1) 



where (f)(xi,X2) = (f){xi + -R, 0^2) is the periodic scalar field, /i ~ [mass] ^"''^''^ is the dimen- 
sional coupling constant and b is the dimensionless parameter of the model. Due to the 
scaling properties of the interaction operator exp(26(/)) + exp(— 260) one can rescale the 
problem to the circle of circumference 2n substituting (6.1) by 



'shG 



27r 

dx2 I dxi 




^2b^ 



+ e 



-2b(j> 



(6.2) 



We are interested in the ground state energy E{R) or, more conveniently, the finite-size 
effective central charge 

6 /? 

Ces{R) = E{R) (6.3) 



TT 



in the ultraviolet limit i? — > 0. Let ^'o[0(a^i)] be the ground state wave functional and define 
again the zero-mode 0o as in eq.(5.6). At — > there is a large region \ogji{R/2TT)'^^'^^ < 
2600 < — log /u(f?/27r)^"'"^^ in the configuration space where one can neglect the interaction 
term in (6.2) and consider 0( free massless field (5.7). In this region the ground 

state wave functional is expected to be a superposition 



*o[0(a:i)] 



(cie2^^'^°+C2e-2^^<^°)|0) 



(6.4) 



of two zero-mode plane waves with some i?-dependent zero-mode momentum P times the 
Fock vacuum |0) of the oscillators in (5.6). The corresponding effective central charge is 
determined at i? — > mainly by P{R) 



Ceff(i?) = l-24p2 + 0(i?2) 



(6.5) 



up to power corrections in R. The momentum P{R) is quantized due to the right and left 
potential walls at 260o ~ ±log//(i?/27r)^+^^ in the action (6.2). Consider say the right 
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wall at 2600 ~ - log//(i?/27r)2+2b . If i? is small enough the second exponent exp(— 26(/)) 
in the potential term of (6.2) is small in this region and does not affect the dynamics which 
is therefore expected to be essentially the same as in LFT. Thus for the reflected wave we 
have 

^o[(j){xi)] ~ {e^'^'f'° + {R/27r)-^'^'^S{P)e-^'^'^°) |0) (6.6) 

where Q is again given by eq.(2.3) and S{P) is the same reflection amplitude (3.18) as 
in LFT (5.12). Factor {R/2tc)~^^^^ appears due to the rescaling of // in (6.2). A similar 
consideration about the left wall reflection leads to the following quantization condition 

(i?/27r)"^^-^^52(P) = 1 (6.7) 

For the ground state momentum this equation reads 

5{P) = TT + 4PQ log(i?/27r) (6.8) 

where we have introduced the reflection phase S{P) 

S{P) = -exp{i5iP)) (6.9) 

Together with (6.5) equation (6.8) determines the most important part of the R ^ 
asymptotic of Cgff (-R). For example, using the regular expansion of the reflection phase in 
the odd powers of P 

S{P) = Si{b)P + S3{b)P^ + S5{b)P^ + ... (6.10) 

one can develop Ceff(-R) systematically in 1/log-R 

^ 247r2 487r%(6) 

cefr(i?) = 1 - + + ... (6.11) 

where we have denoted 

I = 6i (b) - AQ log(i?/27r) (6.12) 

On the other hand the sinh-Gordon model is integrable and its factorizable scattering 
matrix is known [25]. This allows us to compute the same effective central charge Ceff(-R) 
by the thermodynamic Bethe ansatz (TBA) technique [26,27]. In the TEA framework it 
is evaluated as the integral 

Ceff(i?) = j cosh ^ log (l + e-^(^)) de (6.13) 

where e{9) is the solution to the following non-linear integral equation 

if{e-d')\og(l + e-'^^'^J— (6.14) 
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In (6.13) and (6.14) m is the mass of the physical particle in the sinh-Gordon spectrum. 
It is related to the coupling constant /j, as [28] 



ITU 



7(-62) 



"'=r(,^^ir(i + 



40r V2 + 26V y 2 + 262 



^2 N 1 ■2+2b'' 



(6.15) 



The kernel ip{6) in eq.(6.14) contains the information about the sinh-Gordon scattering 
and reads explicitly 

4 sm —r cosh B 

1 + 62 _ 

27r62 

cosh 20 — cos 



1 + 62 

It is straightforward to solve eq.(6.14) numerically. At R small enough eqs.(6.5) and (6.8) 
can be interpreted as a parametric representation of the "experimental" reflection phase 

l-Ces{R) 

24 (6.17) 

^(TBA) ^ ^ ^ log(i?/27r) , 

R being the parameter. According to (6.5) we expect (5'^^^^^(P) to reproduce the Liouville 
phase 5{P) up to exponentially small in 1/P corrections 

S{P) = 5(™)(P) +0 (exp (-2^)) (6.18) 
In particular in the expansion 

5(TBA)(p) ^ 5(TBA)(^)p ^ 5P^)(6)P3 + S^^^^\b)P' + ... (6.19) 

all the coefficients are the same as in eq.(6.10). 

We have solved eq.(6.14) numerically at different values of the parameter 6 and esti- 
mated the few first coefficients ^2k+t\^) (6.19). In Table 1 the numbers are compared 
with the corresponding d2k+i{b) in the expansion (6.10) of the exact Liouville amplitude 
(3.18) 



4 2 ,n^ V2 + 262y V 2 + 262^ , ^ 

5i 6) = -log6 -4giog ^ '—^ - + C 

b 4-y/7r 

16 (6-20) 
'^3(6) = yC(3)(6^ + 6-3) 

64 

S5{b) = -cmb'+b-') 

etc. Here C is the Euler's constant and in 5i (6) the cosmological constant /x is substituted 
in terms of m by means of eq.(6.15). In the numerical calculations we set m = 1. 
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We consider the content of Table 1 as impressive evidence in support of the exact 
Liouville reflection amplitude suggested in sect. 3. The same TBA analysis can be applied 
also for LFT with < 25. Parameter b is complex in this case and the sinh-Gordon 
scattering theory has to be replaced by the staircase model [29]. We hope to say more 
about this interesting relation in a future publication. 

The exact Liouville reflection amplitude can be used in the opposite direction in the 
analysis of the sinh-Gordon (or the staircase) model itself. Subtracting the leading — > 
asymptotic of Ceff(-R) predicted by eqs.(6.5) and (6.8) from the TBA numerical data one 
can separate the power-like corrections in (6.5). Perhaps this would allow us to clarify 
their nature. Work along this line is now in progress. 



7. Conformal bootstrap 

In this section we study numerically the conformal bootstrap equations (2.21) us- 
ing the representation (2.23) of the four-point function in LFT. The structure constants 
C{ai, a2, (Xs) are proposed explicitly in sect. 3. The conformal block JF(Aq,., A, x) which 
also enters eq.(2.23) is not known generally in an analytic form. In refs.[19] the following 
convenient representation has been derived 



J^{Aa,,A,x) = (16?)^'a;^'/^-^«i-^«Kl -a^)'^'/^"^"^"^"^x 
[^3 (g)] +^«^ +^"3 + A.. ^H{XlA\q) 



Here 



and 

is related to x by the equation 



(7.1) 



q = e'""^ (7.3) 



.K{l-x) 



where 

Ki-A = 

'2 Jo ^/t{l-t){l-xt) 
Function H{X'f, A\q) is better parameterized in terms of the variables 

4" - = (I - 

instead of the external dimensions Aq, . . It satisfies the following recurrence relation 

H{Xl A\q) = l + Y^ 9-^^n(A.) ^^^,^ ^^^^ ^ ^^1^^ 
m,n ^ 



X^^9l_A^^^(^-a,) (7.6) 
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Here the sum is over all pairs (m, n) of positive integers and 

g2 (^/6 + n6)^ 

^m,n - ^ ^ 

are the dimensions of degenerate representations of the Virasoro algebra with the central 
charge (2.5). With the notation 

Xm,n = {m/b + nb)/2 (7.9) 
the multipliers Rm,n{Xi) read explicitly 

P X nr,s('^l + A2 — Aj.,s)(Ai — A2 — Aj.,s)(A3 + A4 — Aj.,s)(A3 — A4 — Xr,s) 

(7.10) 

The products in (7.10) are over the following sets of integers (r, s) and {k, I) 

r = — m + 1, — m + 3, . . . , m — 3, m — 1 

(7.11) 

s = —n + 1, —n + 3, . . . , 77. — 3, n — 1 

and 

/c = — m + 1, — m + 2, . . . , m — 1, m 

(7.12) 

/ = — n + 1, — n + 2,...,n — l,n 

while the prime sign near the last product symbol ^ means that the two pairs {k, I) = 
(0,0) and {m,n) are missing. Prom the numerical point of view the expression (7.7) is 
"almost analytic" in the sense that it permits us to compute i/(A|,A|g) very fast and 
precisely. 

In view of eq.(7.1) we find it more convenient to write down the four-point function 
(2.20) as 



and study the crossing properties of the reduced function 

^ J C{ai, a2, Q/2 + iP)C{as, a^, Q/2 - iP) \{mqf' H{\1 Q^/A + P^) 
From eqs.(2.21) it follows that 



(7.13) 



dP 



(7.14) 



^ai,a2,Q:3,Q;4 ("T") ''') — fi'oi,a2,a4,a3 (''' + I5 T + 1) (7.15a) 
9ai,a2,a3,a4{r,f) = \t\^^ '^^i ^'^^ ga^,as,a2,a4{-'^/r, -l/f) (7.156) 
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The first equation (7.15a) is identically satisfied by (7.14) due to the following property of 
the function H{X'^,A\q) 



H{Xl Xl Xl Xl A\q) = H{Xi, Xi, Xi Xi A| - q) 



2 \2 \2 \2 



(7.16) 



which is easily derived from the relation (7.7) and eq.(7.10). Equation (7.15b) still remains 
a non-trivial condition for the structure constants which is believed to contain the whole 
information about the associativity of the operator algebra. 

For numerical calculations we have chosen the correlation function of four puncture 
operators (2.11)^ 



g^j. _ J_ ^ fi'ai,a2,a3,«4 (''"7 ''") 



16 daida2dasda4 
After separating some irrelevant overall factors 



ai=a2=ot3=a4=Q /2 



(7.17) 



9{r,f)- '° 



7r/i7(6 )6 



2\ 1,2-26-^ 



-Q/b 



T — T 

2i 



-QV2 



fir, f) 



(7.18) 



we have to verify the relation 



/(T,f) = /(-l/r,-l/f) 



(7.19) 



for the function 



f{r,f) 



1 f T — T 



2\ 2i 



r{P) {16q)^^H{0,Q'^/4 + P'^\q) dP (7.20) 



Here 



r(P) = 



n^T(^iPyr{-2iP)_ 
TlT^{Q/2 + iP) 

. , 2ttP . , , 
smh — — smh(27roP) exp 




' dt sin^ Pt (1 - e"^* cos2 Pt) 
t sinh(6t) sinh(t/6) 



(7.21) 



We have computed /(r, f) for several values of b and found that eq.(7.19) is satisfied 
up to high numerical accuracy. In figs.l and 2 functions /(t, f) and /(— l/r, — 1/f) are 
compared for purely imaginary r = it at b = 0.8 and b = {1 + i)/V2. This last complex 
value corresponds to cl = 13. 



^ The corresponding values of cti, 0:2, 0:3, 0:4 are well within the domain (2.25) and so 
no "discrete terms" are expected to appear in (2.23) and (7.14). 
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8. Accessory parameters 

The classical Liouville action (2.34) with n singular points is closely related to the 
classic problem of uniformization of Riemann surfaces [30] and in particular to the so- 
called problem of accessory parameters [31]. The last is basically formulated as follows. 
Consider the ordinary linear differential equation 

d^i;{z) + V ( ^. ^ ,„ + i;{z) = (8.1) 

with n regular parabolic singular points Xi. The complex infinity 2; = oo is supposed to be 
a regular point of (8.1) so that the accessory parameters Cj are restricted by the relations 

n 
i=l 



+ 1/4) = (8.2) 

1=1 

n 

5](x2Q + x,/2) = 

The problem is to tune these parameters in such a way that the monodromy group of 
eq.(8.1) is a Fuchsian one. It was proven in ref.[32] that the problem is solved by the 
r/j — > 1/2 version of the Liouville action (2.34). In the case 77^ ^ 1/2 one should be more 
careful since the boundary conditions (2.33) for the Liouville equation (2.30) are slightly 
more complicated 

(^(2;, z) = — 21og |2;|^ + 0(1) at |2;| ^ 00 

(^(2;, 2;) = — log|2; — a^il — 21oglog — a^il + 0(1) at \z — Xi\ ^ 

and more subtractions are needed to regularize the Liouville action (see ref.[32] for more 
details). If S^^^^ {xi, . . . , Xn) is the classical Liouville action of the solution ip{z, z\xi, . . . , x„) 
to (2.29) with the boundary conditions (8.3) then the accessory parameters 

C, = -^S^'^'\xr,...,Xn) (8.4) 

solve the above problem. 

From the LET point of view the action 5'*^'^^'' (xi, . . . , x^) can be considered as the 
leading classical asymptotic 6 — > of the n-point function 



1 d ^ai,...,a„(^l; • • • ; ^n) 



2^ dai . . . dar 



Q;i=...=q;„=Q/2 



e^p(^^^S(-''\x,,...,Xn)^ (8.5) 



of the puncture operators (2.11). In this section we use the representation (2.23) of the 
four-point function to get some information about S^^^^^ {xi, . . . , Xn) and therefore about 
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the accessory parameters in the case n = 4. Note that due to (8.2) in this case there is 
only one independent accessory parameter C{x,x) and eq.(8.1) can be reduced to 

d^ih(z) + + + —, [, ^ [ ih(z) = (8.6) 



where x is defined in eq.(2.19). Eq.(8.5) is reduced to 
1 d'^G (x x) 



16 daida2da3da4 

while 



ai=a2=oi3=a4=Q /2 



exp(^-i5('=')(x,x)) (8.7) 



, dS^''^\x,x) 
C{x,x) = (8.8) 

Let us now consider the classical limit 6 — of eq.(2.23) with four "heavy" operators 
cti = Vi/b- It is convenient to rescale the integration variable P = p/b. The structure 
constants behave as in eq.(3.20) in this limit while the conformal block .^-'(Aq,., A, x) has 
a similar 6 — > asymptotic 

JF(A«„A,a;)~exp(^^/(r7i,p,a;)^ (8.9) 

Function f{r]i,p, x) (which is sometimes called the classical conformal block) is again gener- 
ically unknown in a closed form but can be straightforwardly developed in a power series 
in X 

/(r;„p, x) = iS-S,- 5,) log. + il±ll^Ml±ll^, + 0(.^) (8.10) 

where we have used the notations d = p'^ + 1/4 and Si = r]i{l — rji). At 6 the integral 
in eq.(2.23) is determined by a saddle point p = Ps, i-e., by the minimum of the function 

Sm,m,m,vMx.^) = -^(^'^(r/i, r/2, 1/2 + zp) + ^(^'^(r/a, r/4, 1/2 - zp)- 

f{Vi,P,x) - fiViiP, x) 

where 5'*^'^^'' (?7i, 772, 773) is given explicitly by eq.(3.21). Therefore the classical asymptotic 
of the four-point function is 

Ga^/b,...,a4/bix,x) ~ exp (^-^'5(f^^...,^4(^,^)^ (8.12) 

where 

^riUV2,V3,V4,('^^^^ ~ ^m,V2,V3,V4,{Ps\Xj^) (8.13) 

and Ps is determined by the equation 

d 
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With the exact formula (3.21) it reads exphcitly 



d 

27r - Hog Sr,i,r,2{p) - i^og Srjs,r,Ap) = if {Vi, P, x) + f {rii, p, x)) (8.15) 



where 

^ . ^ r^(l - 2zp)7(?7i + ri2- 1/2 + iphil/2 + m - V2 + ip) 
vi,V2{P) p2(i ^ 2ip)-f{rii + ri2- 1/2 - ip)-f (1/2 + ryi - 772 - ip) 



(8.16) 



The accessory parameter (8.8) corresponds to the special case r]i = r]2 = rjs = r]4 = 1/2 
of eqs.(8.13-16) 

C{x,x)^^f{l/2,ps,x) (8.17) 
At small X one can keep only the leading a; ^ term in the classical block (8.10) so that 

xC{x,x) = {pI-1/4:){1 + 0{x)) (8.18) 
where ps is determined (up to power corrections in x) by the equation 

T(l-2ip)T^(l/2 + ip) 

"'"^^^ + r(i + 2irtr^(i/2-irt = " 

One also can systematically pick up the power corrections in x using the expansion 
/(l/2,p,a;) = (p' - \ogx+ (p' + ]) ^ + 

2 /OQ^2 10 1 \ ^3 



4/ 2 



13p^ + 9 ^ 1 p + ^ 19 ^ 1 p + 



16 32 256(p2 + 1) 7 4 V 24 48 128(p2 + 1) / 8 

(8.20) 



9. Conclusion 

The expression (3.14) for the three-point function, together with the supporting ev- 
idence in sects. 4-8, is the main result of this paper. Although (3.14) is a conjecture we 
find the evidence convincing enough to take it as the starting point in addressing some 
intriguing questions in 2D Quantum Gravity. 

The structure constants (3.14) allow one to construct (in principle) the multipoint 
correlation functions of LFT through the decompositions similar to (2.23). This gives 
access to multipoint correlation fiuictions of Quantum Gravity (coupled to a matter theory) 
at fixed conformal moduli; this is in contrast with the integrated correlation functions 
considered in [11-15]. This moduli dependence of the correlation functions is expected to 
add some insight on the nature of physical states in 2D Quantum Gravity. 

One can also try to analyze the physics of non-minimal CFT (perhaps non-unitary one 
to maintain cm < 1, to begin with) or non-conformal matter QFT coupled to Quantum 
Gravity. 
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The most interesting question is exactly what is it that happens to 2D Quantum 

Gravity when the central charge Cm of matter theory [cm = 26 — c^) exceeds 1. We hope 
that the structure constants (3.14) could be an appropriate vehicle to enter this still rather 
mysterious domain. 

And of course it is important to extend the above analysis to incorporate SUSY. 
We hope to return to this questions in future. 
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Note Added 

After this work was completed we have learned that the conjecture (3.14) already 
exists in the literature. In ref. [34,35] the expression for the Liouville three-point function 
apparantly equivalent to (3.14) is proposed with the motivations very similar to those 
discussed in sect. 3. We are grateful to H.Dorn for bringing the papers [34,35] to our 
attention. 
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Figure Captions 



Fig.l. Reduced four-point function f{T,f) of eq.(7.20) (continuous line) at 6 = 0.8 and real 
t = —IT. It is compared with /(— l/r, — l/r) (points) to verify the crossing symmetry- 
relation (7.19). 

Fig. 2. The same as in fig.l but at 6 = (1 + i)/\/2 corresponding to cl = 13. 
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Table 1. First three coefficients (jC^^^) in the expansion (6.19) obtained by numerical 
analysis the "experimental" reflection phase (6.17) in comparison with the corresponding "exact" 
ones 5*^^^^-' given by eq.(6.20). The uncertainty in the "experimental" numbers is ±1 in the last 
(bracketed) digit. 



l+b2 


^(TBA) 


^(LFT) 


r(TBA) 




r(TBA) 


r(LFT) 


0.1 




-16.61723 




173.3336 




-3225.315 


0.2 


-4.74(4) 


-4.743972 


52.0(5) 


52.08913 


-42(4). 


-425.1404 


0.3 


-0.671(3) 


-0.6715864 


24.64(3) 


24.64884 


-111.(2) 


-111.9785 


0.4 


1.043(1) 


1.042998 


15.26(8) 


15.26739 


-41.4(0) 


-41.39168 


0.5 


1.533(5) 


1.533544 


12.81(9) 


12.82149 


-26.5(6) 


-26.54535 
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